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Abstract — This paper will gives analysis and automatic load 
frequency control of two area power system using Fuzzy Logic 
Controller (FLC) and also with new generation of set of data 
from Genetic Algorithm (GA). Change in frequency is used as 
parental data for getting new generation set of data. Inputs for 
FLC are obtained from the process of Genetic Algorithm 
computational set of data. With improved generation set of data, 
FLC inputs are properly and carefully taken for obtaining 
control vector from defuzzified output of FLC. The tuning of 
FLC is necessary to get an output for both dynamic and static 
performance. The output of self tuned FLC is found reasonably 
well over these conventional controllers. This paper is 
concentrated for developing the algorithm to load frequency 
control and tie line power of two area test system. 


Index Terms — Load Frequency Control, Interconnected 
Power System, Fuzzy Logic Controller and Genetic Algorithm. 
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I. Introduction 

Big importance is giving to control of power generation 
and transmission for obtaining smooth operation of power 
system. The modern power systems with industrial and 
commercial loads need to operate at constant frequency [15] 
with reliable power. Load Frequency Control (LFC) is a very 


important issue in power system operation and control for 
supplying sufficient and reliable electric power with good 
quality. 

Fuzzy [25] base type of controller adds better results in 
system control and also to reducing the steady state error. 
FLC with other techniques like Genetic Algorithm, Neural 
nets and optimization techniques are using for the purpose of 
system manage and control in required area of field. When 
unrestrained case, more oscillation, negative overshoot be 
observed but while comparing to conventional type controller 
P,PI and PID[1], also propose work result gives better 
performances of dynamic responses. 

The interconnected power system is typically divided into 
control areas, with each consisting of one or more power 
utility firms. Sufficient supply for generation of each 
connected area to meet the load demand of its users. Load 
frequency control is one of them. In this work, test system 
with two areas is considered both are used to determine the 
parameters of change in frequency and tie line power 
according to the system dynamics. Both conventional and 
fuzzy are similar in the sense that these two techniques are 
population based heuristic search methods and they approach 
for the optimal solution by latest updates. 


Load 



Figure 1. Block Diagram of test Power System 

For analysis of load frequency control, power system with 
single generator that supplying power to load is considered. 
Block diagram of power system for mathematical model is 
shown in figure 2. 
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Figure 2. Mathematical Modelling of Block Diagram of test 
Power System 
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In order to analyse the load frequency control [4] of power 
system with mathematical modelling, the needful transfer 
functions[15] of model are 
Transfer function of Speed Governor 


1+sT, 





Transfer function of Turbine 

K t 

1 + sT t 

Transfer function of Power System 


1 T sT p 

A. GENERATOR-LOAD MODEL 




Mathematical model of Generator which supply power to load 
AP G -AP d = (— -Af) + (BAf) (4) 


AF(s) = 


APqCsJ - AFfcts) 

B + f - 

r c 


K, 


AF(s) = [AP c (s) - AP D (s)] X ■ _ 

1 T sip 

And we can represent block diagram in figure 3. 


(5) 


(0 



Laplace transform of above equation 
[AP Cl (s) - AF^Cs) - AP tie _JG pL fc) = AF^s) 


(13) 

Block diagram of system is depicted in figure 4. 





AGj.fs] 



Figure 4. Mathematical of Model of Interconnected Power 
System block diagram 

III.DESIGN OF PROPOSED CONTROLLER FOR TEST 
SAMPLE OF POWER SYSTEM 

Simulation of test system can be done with all defined transfer 
functions of block diagram shown in figure 1 and figure 2. 
Original values of data set are change in frequencies AF (s ') 
and AF\s). In many available solutions, GA plays crucial role 
to generate[26] new set of data accordance survival of fitness. 
The evaluation process with development algorithmatic steps 
goes in natural way of related GA theory and concepts. 

To achieve the better set of generation from test data for 
solution has to complete main stages of 

i) Parent Selection and Reproduction 

ii) Crossover 

iii) Mutation 

For getting new generation quietly good and fit for all 
possible test conditions of system environment, the process of 
all stages said above has to be complete with most care and 
systematically. Mathematical representation of data mapping 
and constant for sealing at any instant 


Figure 3. Block diagram of Generator-Load Model 

Steady state response with A P c = 0 and for step load change 

4FW =- ~ 3: - -fe.. - , (7) 


, , r K ' K t K P 

a + sTpJ+ Kl+ S TjCH- S T t ). 


B. INTERCONNECTED POWER SYSTEM 

If losses of tie line[4,15] are neglected, then power transfer 

equation 

IVJ |V 2 | 

A i nil mil I 
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tlE -1 
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5i, 5 2 — Nominal phase angles of voltages and V 1? V 2 at both 
ends. 

Incremental frequency, 


1 d , , 
Af= — — (AS) 
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2ir dt 

A5 - Phase angle deviation 
Incremental tie line power in Area-1, 2 

AP tit _i = 2% T 1Z [J dt- | Af 2 dt p.u. (10) 
AP tit _; = 2 kT 21 If Afj dt - | Af t dt] p.u. (11) 


Net Surplus Power 
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n - Value of data point at instant 
p - Previous data set point 
k max , kmin- User defined constant 
Normal line fit equation 
y = x (k0+k 2 

Data set for new generation was considered A F= AF as 
equalent terms of x and x . Sample test also shown Table 1 
and 2 respectively 

For finite number solutions 

Z 2i - 1 

+ 3 C 15 ) 

i = „ lB n 

- Decreasing sequence of real numbers 
fOO = C - £ [yf- yf~ L ] (16) 

i= 1*2,3._n 

‘C’ can be taken accordance to fitness function scaling. 

Generate set of x with count ‘k th 4 instant 
ms + k f 


= 


fpo m 

Where f = mx+k 

m - data scaling factor with moderation value ‘k’ 
New set of data for updation 


(17) 
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A. GENETIC ALGORITHM FOR REGENERATION OF 
NEW DATA SET 

The sequential steps for searching optimal solution for change 
frequency using GA as shown in figure 5. 



Figure 5. Flow chart of GA Process 

The following algorithm describes the use of GA for 
determining the values of new generation of data sets. 

(i) Select a set of data points from the test system i.e. change 
in frequency of plants. 

(ii) Evaluate the Performance to fit the data in given range of 
input to FLC. 

(iii) With proper selection, off spring crossover and mutation, 
new set of generation of data can obtained. 

(iv) After obtaining the Generation of new individual datasets, 
can be normalized. 

(v) The optimum integral feedback controller gain of case(i) 
and case(ii) at nominal condition. 

The change values of frequencies AF(s') and AF(s') can be 
reproduced with following step by step algorithm.After 
Evaluation process the new data set of generation is tabulated 
in Table 1. 


Data 

No. 

Evaluation 

Function 

Cl 

c 2 

Y, 

Y 2 

f(x) 

m 

5 

1.8 

-3 

2 

-6 

-2 

-0.6 

1.33 

7 

3.2 

-2 

-3 

4 

5 

1.1 

1.12 

16 

2.32 

-6 

2 

4 

3 

1.2 

1.81 

25 

1.93 

-1 

3 

4 

-2 

0.48 

1.21 

31 

2.12 

3 

-1 

2 

6 

0.72 

1.91 

46 

1.62 

1 

4 

3 

-1 

0.61 

1.82 

53 

1.11 

2 

-1 

-1 

1 

0.52 

1.18 

60 

3.02 

6 

5 

1 

-1 

0.92 

1.20 

73 

2.83 

4 

-2 

7 

8 

-1.2 

1.62 

82 

1.06 

4 

1 

-1 

1 

0.29 

1.90 

94 

1.8 

1 

4 

-1 

1 

1.81 

1.02 


Table 1. Sampled data and Evalution of fitness function, f(x). 
Sampled data for reproduction can fit in particular range and 
determine from evaluation function using equation 15,16. 
Two cross over point[22] used to get new generation. New 
generation set of data is shown in figures 6, 7 for 100 and 200 
iterations 


New Generation of set of DATA-100 iterations 



Figure 6. New generation set of data-100 iterations 

New Generation of set of DATA-200 iterations 



Figure 7. New generation set of data-200 iterations 
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B. FUZZY LOGIC CONTROLLER WITH GA 
REPRODUCED INPUTS 

The main components of fuzzy control system are the 
following blocks: 

- Fuzzification, 

- Rule base, 

- Inference engine and 

- Defuzzification. 

Merits of fuzzy logic controllers are mainly solves the non 
linear problems with vague inputs, not needed an accurate 
mathematical model, handling nonlinearity. For controlling 
such a complicated system, fuzzy logic controllers gives very 
good results like this application. Inputs for fuzzy logic 
controller are used from new generation set of data as 
obtained from GA process. 



Figure 8. Fuzzy Logic Controller with New generation set of 

data as input 

Inputs to a fuzzy logic controller are usually an error ‘e’ and a 
change of error ‘Ae\ Here the two inputs for fuzzy logic 
controller are AF(s) and AF'iV). The values that indicate 
input variables change of frequency and derivative of change 
in frequency. Linguistic terms used for the membership 
functions are such that are terms NL(Negative Large), 
NM(Negative Medium), NS(Negative Small), ZE(Zero), 
PS (Positive Small), PM(Positive Medium) and PL(Positive 
Large). 

Fuzzification block converts crisp inputs to one or many 
membership grade values; where the change in frequency and 
derivative of change in frequecy accordance with variation of 
tie-line power, are described by membership functions given 
in Figure 9 and Figure 10. Membership function of output of 
control vector is shown in Figure 11. Afterwards, it is possible 
to apply descriptive rules of reasoning Power searching was 
POSITIVE and LARGE and the last change of desired change 
in frequency was PL then keeps tracking the derivative of 
change in frequency in the same PL direction with PL 
increment. Rules like this are involved in block “rules table”, 
and they are given in Table 2. Finally, the fuzzy set of output 
reference change in frequency is back “deffuzified” to convert 
it to the actual value. 

That means the output values such are PL, PM, ZE, 
SMALL are translated to numbers which indicates a 
measurable (but normalized) value of the change in 
frequency. It could be also noticed the output of controller is 
added by some amount of in order to avoid local minima in 
characteristics due to the changes in frequency value. That 
means the output values such are PB, PM, PS are translated to 
numbers which indicates a measurable (but normalized) value 
of the changes in frequency. Min-max implication method has 
been used for getting output from all rules of decision table. 




Figure 10. 2 nd Input membership functions, AF(s’) 



Figure 11. Output Membership Function, u 
Fuzzy control is based on a logical system called fuzzy 
logic. It is much close in spirit to human thinking and natural 
language than classical logical systems. The LFC as has been 
reported in several papers is to maintain balance between 
production and consumption of electrical power. Due to the 
complexity and multi-variable nature of power systems, a 
conventional control method has not provided satisfactory 
solutions. The fuzzy logic control has tried to handle the 
robustness, reliability and nonlinearities associated with 
power system controls. 

The block diagram of figure 2 shows the single area fuzzy 
load frequency controller. It uses two input membership 
variables (frequency deviation AF and derivative of change in 
frequency deviation AF'fsj). Fig. 9, 10 and 11 shows the 
membership functions for the AF , AF(sj and the control 
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output u. The input signals are first expressed in some 
linguistic variables using fuzzy set notations such as 
NL(Negative Large), NM(Negative Medium), NS(Negative 
Small), ZE(Zero), PS (Positive Small), PM(Positive Medium) 
and PL(Positive Large). 
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AFts) 

NL 

NM 

NS 

ZE 

PS 

PM 

PL 

NL 

ZE 

PS 
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PL 

PL 

PL 

PL 

NM 

NS 
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NL 
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PS 

PM 

PM 

NL 

NL 

NM 

NM 

NS 

ZE 

PS 

PL 

NL 

NL 

NL 

NL 

NM 

NS 

ZE 


Table 2: Decision table of 7x7 rule base 

Table 2 shows a set of decision rules [9], also expressed in 
linguistic variables relating input signals to the control signal. 

The rules interpreted as follows: 

1. If (AF(s) is NL) and (AFi's) is NL) then (Control Vector, u 
is ZE) 

2. If (A F ( s j is NL) and (AF(s) is NM) then (Control Vector, u 
is PS) 

3. If (AF(s) is NL) and (AF(sj is NS) then (Control Vector, u 
is PM) 

4. If (AF(sj is NL) and (AFi'sj is ZE) then (Control Vector, u 
is PL) 

5. If (AF (sj is NL) and (AF(Vi is PS) then (Control Vector, u 
is PL) 

6. If (AF(sj is NL) and (AF(Vi is PM) then (Control Vector, u 
is PL) 

7. If (AF(sj is NL) and (AF(Vi is PL) then (Control Vector, u 
is PL) 

8. If (AF(sj is NM) and (AF"(sj is NL) then (Control Vector, u 
is NS) 

9. If (AF(sj is NM) and (AFl(s) is NM) then (Control Vector, 
u is ZE) 

10. If (AF(sj is NM) and (AF'(s) is NS) then (Control Vector, 
u is PS) 

11. If (AF(s) is NM) and (AF(Vi is ZE) then (Control Vector, 
u is PM) 

12. If (AF(sj is NM) and (AF’(sj is PS) then (Control Vector, 
u is PM) 

13. If (AF(s) is NM) and (A¥\s]\ is PM) then (Control Vector, 
u is PL) 

14. If (AF(s) is NM) and (AF(Vi is PL) then (Control Vector, 
u is PL) 

15. If (AF(sj is NS) and (AF(s) is NL) then (Control Vector, 
u is NM) 

16. If (AF(sj is NS) and (AFfVi is NM) then (Control Vector, 
u is NS) 

17. If (AF(sj is NS) and (AF(Vi is NS then (Control Vector, u 
is ZE) 

18. If (AF(s) is NS) and (AFi's) is ZE) then (Control Vector, u 
is PS) 


19. If (AF(s) is NS) and (API's) is PS) then (Control Vector, u 
is PS) 

20. If (AF(sj is NS) and (AFi'sj is PM) then (Control Vector, 
u is PM) 

21. If (AF(sj is NS) and (AF(s) is PL) then (Control Vector, u 
is PL) 

22. If (AF(s) is ZE) and (AFl'sj is NL) then (Control Vector, u 
is NM) 

23. If (AF(sj is ZE) and (AF(Vi is NM) then (Control Vector, 
u is NM) 

24. If (AF(s) is ZE) and (AFi'sj is NS) then (Control Vector, u 
is NS) 


Similarly remaining all rows of decision table will be 
written as TF’ and ‘THEN’ propositions. Propositions (IF, 
THEN) with all 49 rules are used to develop programming 
code in MATLAB. After developing rule base in matlab, the 
following figure 12 show all rules [6] in MATLAB window. 
Ligure 12 interpretation of some sample of rules with input 1 
is 0.349 and input 2 is 0.50 and also evaluation of these two 
inputs to obtain control output vector. The value of control 
vector, u is determined as 0.266 from simulation of two inputs 




Rule Editor LoadFrequencyControl 





File Edit View Options 


1. II (Change in Frequency is NL) and (Derivative of Change in Frequency is NL) then (Control Vector is ZE) > 

2. If (Change in Frequency is NL) and (Derivative of Change in Frequency is NM) then (Control Vector is PS) 

3. If (Change in Frequency is NL) and (Derivative of Change in Frequency is NS) then (Control Vector is PM) - 

4. If (Change in Frequency is NL) and (Derivative of Change in Frequency is ZE) then (Control Vector is PL) 

5. If (Change in Frequency is NL) and (Derivative of Change in Frequency is PS) then (Control Vector is PL) 

6. If (Change in Frequency is NL) and (Derivative of Change in Frequency is PM) then (Control Vector is PL) 

7. If (Change in Frequency is NL) and (Derivative of Change in Frequency is PL) then (Control Vector is PL) 

8. If (Change in Frequency is NM) and (Derivative of Change in Frequency is NL) then (Control Vector is NS) 

9. If (Change in Frequency is NM) and (Derivative of Change in Frequency is NM) then (Control Vector is ZE) 

10. If (Chanqe in Frequency is NM) and (Derivative of Chanqe in Frequency is NS) then (ControiVector is Pi 

<1 m I ► 


If 


Change in 


and 

Derivative of 


Then 

Control Vector is 




> 

NM 


NS 

^ | 

ZE 


PS 


PM 





NS 

> 



PS 


PM 


PL 


none 

▼ 


not 


Delete rule 


Add rule Change rule 


d-d 


FIS Name: LoadFrequencyControl 



Ligure 12. Rule base in Matlab window 
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After developing rule base in matlab, the following figure 13 
gives show the how rules will excute the two inputs to all 

propositions. 

IV.SIMULATION OF TEST SYSTEM AND RESULTS 


Simulation is carried for two area test power system. Each 
area is represented by governor, turbine and plant; also two 
areas are combined with tie line. MATLAB simulink is used 
for development of test system and shown in figure 14. 



Figure 14. Matlab Simulink model of two area power system. 


Change in frequency of area-1, AFi and change in frequency 
of area-2, AF 2 are taken to determine change in frequency AF. 
To test power system transfer function blocks are used in 
simulation model. Obtained transient response of AF 1? AF 2 are 

shown figures 15, 16 for simulated test system. Output 
response of tie line power is shown in figure 17. 


0.39 

0.38 

0.37 

0.36 

0.35 

0.34 

0.33 

0.32 

1 23456789 10 

Time in seconds 

Figure 15. Frequency deviation in area-1, AFy 


Variation of Change in frequenccy for areal 



Variation of Change in frequenccy for area2 



Figure 16. Frequency deviation in area-2, AF : 


x 10' 


Variation of Change in tie line power 



Figure 17. Change in tie line power, AP tie 
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Simulation carried out for response of load variation (1-2%) 

and shown in figures 18 to 23. 

Variation of Change in frequenccy for areal 



Figure 18. Response of AF { for 1% load variation 


Variation of Change in frequenccy for area2 



Figure 19. Response of AF 2 for 1% load variation 


x 10 


-4 


Variation of Change in tie line power 



Figure 20. Response of AP tie for 1% load variation 


Variation of Change in frequency in areal 



Figure 21. Response of AF! for 2% load variation 

Variation of Change in frequency in area2 


[ L L [ 


. WITHOUT CONTROLLER 

-PROPOSED CONTROLLER 



0.4 



Q 2 L _[_[_[_[_[_[_[_[___[ 

1 23456789 10 

Time in seconds 


Figure 22. Response of AF 2 for 2% load variation 


Variation of Change in tie line power 



Figure 23. Response of AP tie for 2% load variation 

Comparative study of output response (Change in 
frequenciesAF 1? AF 2 and change in tie line power AP Tie ) is 
shown in table 3 and table 4. Also output response (Change in 
frequenciesAF l9 AF 2 and change in tie line power AP Tie ) for 1 
to 2% load variation also shown in table 3 and is shown in 
table 4. 
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AF-Areal 
(Sec) 

AF-Area2 

(Sec) 

APTie 

(Sec) 

Without 

Controller 

Without 

load 

variation 

10 

9 

9 

1% load 
variation 

13 

8 

10 

2% load 
variation 

9 

11 

11 

Proposed 

Controller 

Without 

load 

variation 

9 

7.5 

8 

1% load 
variation 

8 

7 

8 

2% load 
variation 

7.2 

10 

9 


Table 3. Results Comparison-Settling time 




AF-Areal 
(p.u) 

AF-Area2 

(p.u.) 

APxie 

(p.u.) 

Without 

Controller 

Without 

load 

variation 

0.385 

0.388 

0.0009 

1% load 
variation 

0.95 

0.35 

-1.2 

2% load 
variation 

0.055 

0.0875 

3.2xl0' 4 

Proposed 

Controller 

Without 

load 

variation 

0.371 

0.371 

0.0006 

5 

1% load 
variation 

0.81 

0.30 

-0.6 

2% load 
variation 

0.05 

0.0825 

2.1xl0' 4 


Table 4. Results Comparison - Overshoot 


V. CONCLUSION 

This paper has presented a control strategy for improving 
performance and stability of two area power system based on FLC 
and new generation set of data calculated from Genetic Algorithm 
process. In proposed work the inputs to FLC are calculated from 
fitness function of genetic algorithm. Frequency deviation in each 
control area and tie line power deviation of test system has been 
considerably increased with proposed method. The simulation 
studies show that proposed method can able to reduce the 1-2% load 
variations of two area test system. 

APPENDIX 

Base and rated values of Generator are as follows 


Base frequency, f base - 50Hz 

Base voltage, V base - 550V 

Base MVA, S base - 5MW 

Base angular frequency, co base - 2 tt rad /sec 

Stator resistance - 0.021 p.u. 

Rotor resistance - 0.015 p.u. 

Rated Voltage - 0.39 kV 

Rated Power - 3MW 

Stator/Rotor ratio - 0.36 

Angular Moment of Ineteria, J -1.5 p.u. 
Mechanical damping - 0.01 

Stator leakage inductance - 0.54 


Rotor leakage inductance - 0.23 
Mutual inductance - 4.2 

Typical values transfer function parameters of Control 
area-1are 

T gl =0.02, T tl =0.09, H 1= 0.132, R 1= 2.1, D 1= 0.0069, T 12 =0.52. 
K gl =5, K tl = 0.5, K pl =20. 

Typical values transfer function parameters of Control area-2 
are 

T gl =0.013, T tl =0.05, H 1= 0.52, R 1= 1.0, D 1= 0.021, T 12 =0.31. 
K gl =4, K tl = 1.8, K pl =10. 
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